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The production of the charmonium states in hadronic experiments is considered 
t^- ' in NLO at the strong coupling constant. It is shown, that such an approach solves 

o • 

some significant problems, arising when only leading order processes are considered. 
In particular, in such a consideration distributions over the transverse momentum of 



oo 



o 
m 
o 

o 



the final charmonium can be obtained. There appears also a natural explanation to 
the existence of Xci-meson in final state, that is observed experimentally and cannot 
be produced in leading order processes. 



INTRODUCTION 



It is well known, that the specific property of the charm production processes is high 
Qh! sensitivity to the gluonic contents of the interacting particles. This property is used not 

o . rj 

only in the determination of gluon distribution in the proton G(x) but also, in the case 
of polarized beams, for the determination of the gluon polarization fraction AG/G [2]. That 
is why new experiment, labeled as SPASCHARM, on hadronic production of charmonia was 
recently proposed. The goal of this experiment is detailed analysis of charmonia production 
in polarized proton collisions at the 70 GeV energy. 

The important part of the Program is the detailed analysis of all possible mechanisms 
of charmonia production. Such an analysis is especially important since at low energies the 
contributions of gluon-gluon, quark-gluon and quark-antiquark subprocesses are comparable. 
For example, if the energy of the proton beam is equal to 40 GeV, the ratio of ip production 
cross sections in pp and pp collisions equals cr(pp)/a(pp) ~ 6, so it is clear that the quark- 
antiquark annihilation plays an important role. As the energy rises this ratio tends to unity 



*Electronic address: Likhoded@ihep.ru 
'''Electronic address: Alexey.Luchinsky@ihep.ru 



2 



and the gluon-gluon process dominates. 

Another problem is that the direct production of -^-meson is suppressed in comparison 
with the production of the intermediate P-wave states Xco,i,2 with the subsequent decay 
Xc — * J/'tpl- m the experiments {3] this fact is confirmed with one important remark. The 
experimentally observed cross sections of Xc2 and Xci are comparable (xecr m eson can hardly 
be observed due to its small radiative width), while the well known Landau- Yang theorem 
forbids the formation of the axial meson from two massless gluons. Thus, the experimental 
observation of this state shows that virtual gluons give additional contribution that is not 
accounted in the usual description of gluon distributions. There are also large uncertainties 
in the available parameterizations of these distributions in the region of small virtualities 
and the values of the gluon momentum fraction x g G [0.1; 0.5] [l|]. Due to these uncertainties 
large errors in the prediction of the ip and Xc production cross sections appear. One more 
difficulty is that the distributions G(x) are integrated over the transverse momentum. As a 
result, such method does not allow one to obtain the distributions of Xco~ and Xc2-mesons 
over pt- 

Initially this problems were solved by the introduction of color-octet (CO) components 
of the quarkonia, that arise naturally in the non-relativistic quantum chromodynamics 
(NRQCD), where the expansion over the relative velocity of quarks in meson is performed. 
In this model it is assumed, that final meson is formed from heavy quark pair in color-octet 
state, that subsequently transforms into a physically observed colorless meson. In the frame- 
work of NRQCD the probabilities of these transitions are described by the matrix elements 
of four-fermion operators, that are determined from the experimental distributions over the 
transverse momentum of final charmonium. In the works {4, 5] it was shown, that octet 
components give satisfactory agreement with the experimental data on TEVATRON at the 
energy yfs = 1.8 TeV. We would like to stress, however, that this explanation will not work 
for charmonium production at lower energies. The reason is that the distributions caused 
by octet components decreases slowly with the rise of the transverse momentum, but the 
probability to find such a component in the meson is small, compared with the singlet case. 
As a result, in the large transverse momentum region the contribution of octet components 
can be significant, but for small energies and transverse momenta it is suppressed. 

Recently another way to solve this problem was proposed, where the so called non- 
integrated over the transverse momentum distribution functions G(x, kx) are used (fey- 
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factorization) . In this case both mentioned above problems are solved simultaneously 

. The transverse momentum of the produced in gluon fusion Xco,2 mesons is explained by the 
transverse momenta of the initial partons. Axial charmonium meson can also be produced 
in gluon fusion, since in the framework of kx factorizations gluons have non-zero virtuality 
of the order k\. There is a number of works, that explain the experimental distributions on 
TEVATRON with the help of these functions (see, for example E], H])- According to these 
works, there is no need to introduce CO components to reproduce the experimental data 
on P-wave charmonium px distributions. Thus, in /cy-factorization approach color-singlet 
components give the main contribution. 

Unfortunately, the method, used in the modeling of the unintegrated distribution func- 
tions G(x, kx), is based on the summation of large log(l/x), so it is not applicable for low 
energies, where the gluon momentum fractions are in the range 0.1 < x g < 0.5. For this 
reason we are forced to use the following approximation in our calculations. We start from 
the collinear gluon distributions with well known collinear distribution functions. Further 
we consider the charmonia production at next to leading (NLO) order in the strong coupling 
constant a s . Such a trick enables us to obtain the distributions over px for all charmonium 
states. For Xco an d Xc2 production we observe a collinear singularity at px = 0. To avoid this 
singularity we introduce a cutoff on the transverse momentum, and the value of the cutoff 
parameter is determined from the inverse geometrical size of the charmonium meson. For 
directly produced ip and Xa such a singularity is absent and we use the whole integration 
region for px- 

In the next section we will briefly describe the formalism used in the subsequent paper. 
Section 3 is devoted to the consideration of different modes of charmonia production and 
analytical expressions for corresponding partonic cross sections are presented. In the fourth 
section we determine the cross sections of the hadronic processes and present numerical 
results. In the last two sections we give our estimates for spin asymmetries and briefly 
discuss the results of the paper. 
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2. CHARMONIA 

One of the characteristic features of charmonium states is the smallness of the relative 
quark velocity v: 

v 2 ~ [a s {m c v)f « 0.2. 

Due to the existence of this small parameter the processes of charmonia production actually 
evolves in two almost independent steps: 

• hard process, in which the quark- ant iquark pair is created, 

• subsequent hadronization of this pair into an experimentally observable meson. 

Second characteristic feature of charmonia is the smallness of the strong coupling constant 
at the hard step of their production. So it is possible to use the perturbation theory while 
calculating the hard part of the amplitude. 

The second stage of charmonia production, i.e. the hadronization of the quark- ant iquark 
pair into a final meson with the momentum p and mass M x \ is described by the follow- 



ing simple procedure First we use the common QCD to write the amplitude of the 



production of on-shell quark and antiquark with momenta p + q/2 and p — q/2 respectively: 

M h ^ d = u{ P - + q)Mv{ P --q). 

Here u and v are the wave functions of quark and antiquark, p is the momentum of the final 
meson (p 2 = M 2 ), and q is the relative momentum of quark-antiquark pair (pq = 0). Then 
we project this amplitude to the color singlet state with the value of the total spin S = 1: 



y 

where i and j are spinor indices of quark and antiquark, m = M/2 is c-quark mass, and 
is spin polarization vector of the pair. The resulting amplitude is then expanded into a 
series in relative momentum q: 

M = e% (Mn + M^q" + ...). (1) 



J ) In what follows we will neglect the difference between tjj and Xc masses 
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The amplitude of the S-wave charmonium meson production can be derived from the first 
term of this expansion: 

where e M is the polarization vector of if)- meson, and R^(0) is the radial part of its wave 
function at the origin. Amplitudes of P-wave charmonia production are derived from the 
second term of the expansion ([!]): 



1 1^(0)1 
V3 4?rm 



\R' 



2V2m 4vrm ^ 



4:71771 
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Here R' x (0) is the derivative of the radial part of \- meson wave function at the origin, e a x 
and e£ CT are polarization vector and tensor of XcV an d Xc2" m esons, and the tensor J pu is 

~ M 2 

The value of radial part of charmonium wave function can be determined by solving the 
Schrodinger equation, or from the experimental values of the decay widths of these mesons. 
In our paper we will use the last way. ^-meson wave function is determined from its leptonic 
decay width: 

The derivative of the radial part of Xco,2-meson wave function can be determined form the 
total widths of these mesons, that are approximately equal to widths of the decays into a 
pair of massless gluons: 

r (Xo) = 96a s 2 ^^, (2) 

x 

128 2 |P'(0)| 2 ^ 

x. 



The ratio of these width 



r(xco) 15 
r(x c2 ) 4 • 
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Figure 1: gg -> XcJ 



agrees well with the experimental data. 

It should be noted, that in the described above procedure we have neglected the relative 



motion of quar ks in t 



was shown 
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re charmonium and use the Bethe-Heitler approximation. Recently it 
Q, 0]N that relative motion of quarks results to significant increase 



of exclusive cross section. It is still not clear, whether it is valid to neglect this motion in 
the calculation of inclusive cross sections. 

The other approximation used in our paper is the CS-approximation 



17 



18 



19 



20] 



In other words we have considered only CS contributions and neglected CO components 
of charmonium mesons. As it was mentioned in the introduction, we think that such an 
approach can be used, because the experimental data on the charmonium production on 
TEVATRON is reproduced in the framework of this model 



2.1. 



gg -> Xc 



In the following we will need the expressions for effective vertex, that describes the tran- 
sition of gluon pair into a Xcj state. The corresponding diagrams are shown on fig. [TJ Using 
the presented above technique we have derived analytical expressions for these vertices. In 
the case of scalar meson, for example, we have 

Ay^a s \R'(0)\ 



A {0) 



V3M3( qi q 2 ) 2 



(3M 2 



it 



y 2n0 



9 



a/3 



(3M 4 -4( gi 2 + g 2 2 )M 2 + ( gi 2 -g 2 2 ) 2 ) . 



(4) 
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The expressions for axial and tensor mesons are rather tedious, so we will not present them 
here. 

Squaring these amplitudes and averaging over the gluon polarizations we have 



where the explicit form of the functions B^ J \q 2 , q 2 ) is 



8Tra 2 s \R' x (0)\< 
3M3( gi g 2 ) 4 

+ 28g 2 ^M 4 - 8q^M 2 - 8q%M 2 + 8q\q\M 2 + 8q\q\M' + q\ + q\ 
- 4g 2 g 2 6 + 6g 4 g 4 -4g^), (5) 



£ (0) (?U 2 ) = 8 3M3^ ( T { 9M8 " 24g ' M6 ~ 2AqlM " + 22qtM4 + 22qtMA + 



B^\ql ql) = ^^^ {gi 2 ^ + ^ - ^ - ^ " U ^ M " ~ 

- q\M 2 - q\M 2 + 17g 2 g 4 M 2 + 17g 4 g 2 M 2 + q\ + g 2 8 - 2g 4 g 4 } (6) 

^ (2) te 2 ^2) = 16 W3^ ( T | 6M8 - - ^ + ^ + ^ + 

?>M A {q x q 2 y I 

+ 34g 2 g 2 M 4 + q\M 2 + g 2 6 M 2 - g 2 g 4 M 2 - q\q\M 2 + qf + q 8 2 - 

~ 4g 2 g 2 6 + 6g 4 g 4 -4^ 2 }. (7) 

It can be seen that in the case q\ = q 2 = we have B^ = 0, i.e. an axial x-meson can not 
be produced from the fusion of two massless gluons, exactly as the Landau- Yang theorem 
states. This is not true for scalar and tensor mesons, so the decays of these mesons into two 
gluons is possible. The widths of these decays can be expressed through the functions B^ 
by the relation 

that directly leads to eq. (TSJ) and (jHJ). One other interesting property of B^ J '{q 2 , q 2 ) functions 
is that in the case of axial meson we observe a strong dependence on the virtualities of the 
gluons, while for scalar and tensor mesons this dependence is not so crucial (see fig. [2]). 



3. PARTONIC SUBPROCESSES 



In this section we will consider the partonic subprocesses that give the contributions to 
production of charmonium Q in the hadronic experiments. At the leading order (LO) of the 
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Figure 2: B^ J \ql = -lGeV 2 ,q%) 

perturbation theory (i.e. ~ a 2 s ) the only possible reactions are the processes gg — > Q with 
Q — Xco or Xc2- Vector charmonium production is forbidden by the charge parity conser- 
vation, and Q = Xci cannot be produced because of the Landau- Yang theorem. However, 
from the experimental data we know, that the cross sections of x& an d Xc2 production are 
comparable j^. The only way to solve this problem is to consider next to leading order 
processes (NLO) in the strong coupling constant. For this reason we will deal with the 
following partonic processes: 

a{h)b{k 2 ) - Q(p)c(k 3 ), (8) 

where in the parenthesis the momenta of the particles are shown. In the case of proton- 
proton interaction we will consider the reactions 

99 -> Qg, 19 -> Qq, 

where q is the valence u or (i-quark. In the case of pp-annihilation the quark- ant iquark 
annihilation is also possible 



qq -> Qg. 
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Figure 3: gg -> 



Consideration of NLO in strong interaction constant solves also the other mentioned in 
the introduction problem. In the widely used collinear approximation the integrated over 
the transverse momentum kr distribution functions are used. As a result, it turns out that 
initial partons and, hence, final charmonium, do not possess any transverse momentum. 

In the rest of this section we will consider in detail the mentioned above partonic reactions 
and give the analytical expressions for the corresponding cross sections. 



Because of the charge parity conservation color-singlet ^-state cannot be produced in the 
reaction gg — > x/j, so at LO it can be produced only with a gluon in the final state. One of 
the diagrams describing this process is shown in figure [3j others can be obtain from it by 
gluon permutations. Using the formalism described in the previous section we have obtained 
the following expression for the partonic cross section: 



3.1. gg -» Qg 



3.1.1. Q = V 




where the definitions for Mandelstam variables s, i and u are 



s = 



{h + hf, i={h-h) 2 



u = {k\ — k 2 ) 2 . 
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(b) 



(c) 



Figure 4: gg -> XcJff 



These variables are not independent, the following relation between them exists: 



s + t + u = M' 



Up to notations our result coincides with the result of work 221 ] . 



3.1.2. Q = Xco,2 



In the case of Xco,2" m eson productions diagrams shown in figure [3] give contribution to 
partonic cross section (obviously, if) in these diagrams should be replaced by Xc)- It is also 
necessary to include diagram shown in figure HI The partonic cross sections of these processes 



tu s 3 (s + *) 4 ( s + m) 4 (* + m) 4 M 3 L 

- 6(* + u) 3 (9P + 14m* + 9m 2 ) M 18 + 

+ (* + m) 2 (153* 4 + 492m* 3 + 695m 2 * 2 + 492m 3 * + 153m 4 )M 16 - 

- 2(* + m) 3 (135* 4 + 393m* 3 + 545m 2 * 2 + 393m 3 * + 135 m 4 )M 14 + 
+ 2(162* 8 + 1065m* 7 + 3208m 2 * 6 + 5852m 3 * 5 + 7096 m 4 * 4 + 

+ 5852m 5 * 3 + 3208m 6 * 2 + 1065m 7 * + 162m 8 )M 12 - 

- 2(135 * 9 + 966m* 8 + 3215m 2 * 7 + 6627m 3 * 6 + 9351m 4 * 5 + 

+ 9351m 5 * 4 + 6627m 6 * 3 + 3215m 7 * 2 + 966m 8 * + 135m 9 )M 10 + 
+ (153* 10 + 1170 m* 9 + 4249m 2 * 8 + 9722m 3 * 7 + 15548m 4 * 6 + 
+ 18124m 5 * 5 + 15548m 6 * 4 + 9722m 7 * 3 + 4249m 8 * 2 + 1170m 9 * 
+ 153m 10 ) M 8 - 2(27? 1 + 222m* 10 + 885m 2 * 9 + 2237m 3 * 8 + 
+ 4001m 4 * 7 + 5308m 5 * 6 + 5308m 6 * 5 + 4001m 7 * 4 + 2237m 8 * 3 + 
+ 885m 9 * 2 + 222m 10 * + 27w n )M 6 + (9* 12 + 90m* 11 + 416m 2 * 10 
+ 1190m 3 * 9 + 2394m 4 * 8 + 3582m 5 * 7 + 4090m 6 * 6 + 
+ 3582m 7 * 5 + 2394m 8 * 4 + 1190m 9 * 3 + 416m 10 * 2 + 90m 11 * + 
+ 9m 12 )M 4 - 2*m(* 2 + ui + m 2 ) 2 (3* 7 + 15m* 6 + 37m 2 * 5 + 
+ 55m 3 * 4 + 55m 4 * 3 + 37m 5 * 2 + 15m 6 * + 3m 7 ) M 2 + 
+ * 2 m 2 (* + m) 2 (* 2 + ui + m 2 ) 4 ] , 
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da(gg^ X c2g) 1 4a>|i^(0)| 



12(t + m) 4 M' 
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dt tMS 3 (s + t) 4 ( s + u) 4 (t + w) 4 M 3 L 

- 24(t + m) 3 (3? + hut + 3m 2 ) M 18 + 

+ (i + m) 2 (204? + 651m? + 880m 2 ? + 65lM 3 t + 204m 4 )M 16 + 
+ (-360? - 1995m? - 4949m 2 ? - 7428m 3 ? - 7428m 4 ? - 

- 4949m 5 ? - 1995u 6 * - 360m 7 ) M 14 + 

+ (432? + 2526u? + 6652m 2 ? + 10877u 3 ? + 12640m 4 ? + 10877m 5 ? + 
+ 6652m 6 ? + 2526m 7 £ + 432m 8 )M 12 + 

+ (-360? - 2274m? - 6290m 2 ? - 10647m 3 ? - 13185m 4 ? - 13185m 5 ? - 

- 10647m 6 ? - 6290m 7 ? - 2274u 8 t - 360m 9 )M 10 + 

+ (204?° + 1455m? + 4328m 2 ? + 7504m 3 ? + 9232 m 4 ? + 9614m 5 ? + 

+ 9232u 6 ? + 7504m 7 ? + 4328m 8 ? + 1455u 9 t + 204u 10 )M 8 + 

+ (-72? 1 - 615u?° - 2085u 2 ? - 3878m 3 ? - 4748u 4 ? - 4678u 5 ? - 

- 4678u 6 ? - 4748u 7 ? - 3878m 8 ? - 2085u 9 ? - 615m 10 £ - 72u n )M 6 + 
+ (12? 2 + 144m? 1 + 616 m 2 ?° + 1345m 3 ? + 1824m 4 ? + 1806m 5 ? + 

+ 1688m 6 ? + 1806m 7 ? + 1824m 8 ? + 1345m 9 ? + 616m 10 ? + 144m 11 t + 
+ 12m 12 )M 4 - £m(? + ui + u 2 ) 2 (12? + 60u? + 
+ 91m 2 ? + 49m 3 ? + 49m 4 ? + 91u 5 ? + 60uH + 

+ 12m 7 )M 2 + 2?m 2 (£ + u) 2 (? + m£ + m 2 ) 4 ]. (10) 

It should be noted that these cross sections diverge in the regions t -> and i -> M 2 -s (that 
is m — > 0). Such divergency is caused by the fact that in these regions the propagators of the 
virtual gluons shown in fig. H] become zero. As a result the cross sections of Xco,2"P r oduction 
at t — > can be written in the factorized form: 

^ (!M ;^ g) - WflTBMftO)- (11) 
at t-»o t 

Here the first term 

f 384rra s (g 2 - gM 2 + M 4 ) 
9 WW ~ M 4 s(s-M 2 ) {lZ) 

does not depend on the spin of the final charmonium and describes the splitting of the initial 

gluon into gg* pair, while the second term describes the production of Xcj meson [functions 
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B(J) were presented earlier, formulas (jSJ) — (CO)]- 

It is easy to understand the physical nature of this divergency. In the regions t — > and 
u — > the transverse momentum of the final charmonium 

(13) 

tends to zero and it is impossible to distinguish the reaction ([8]) from reaction of Xc produc- 
tion in the gluon pair fusion, that is the process gg — > Xco,2- To avoid this divergency we 
will impose the following restriction on p?: 

p T > A 

and use the inverse geometrical size of the charmonium as the cutoff parameter A: 

1 

A ~ ~ m c v . 

Rcc 

It is interesting to mention, that the ratio of this divergent cross sections tends in the 
region of small pt to a finite value, that is equal to 

lim f d °(99 -> Xc2g) j dajgg -> Xcog) 1 = 2J X2 + 1 r(x e2 ) = 4 
pt-o\ dp T I dp T J 2J X0 + ir(x c0 ) 3' 

This relation follows immediately from the eq. ffTTj) . since in the range of small pt we have 
the production of Xc-mesons from two almost massless gluons. That is why the partonic 
cross sections are, up to spin factor, proportional to the widths of the decays of these mesons 
into a massless gluon pair. 
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3.1.3. Q = X 



cl 



If Q = Xci all shown in figures [3l H] diagrams give contributions to the partonic cross 
section, and the expression for this cross section is 

dv{gg -> xag) 



12a 3 vr|i4| 



dt 



(t + u) 2 (t 2 + u 2 )M 



ii 



s 2 (s + t) 4 (s + u) 4 (t + w) 4 M 3 

- 4(2i 5 + 5wt 4 + 6w 2 t 3 + 6u 3 i 2 + 5uH + 2u 5 )M 12 + 

+ (26t 6 + 80wt 5 + 115m 2 t 4 + 120uH 3 + 115w 4 P + 80w 5 t + 26w 6 )M 10 - 

- 2(22t 7 + 82ui 6 + Ulu 2 P + 164m 3 £ 4 + 164m 4 F + 141w 5 t 2 + 82uH + 
+ 22m 7 ) M 8 + (Alt 8 + ISAui 7 + 378w 2 £ 6 + 510w 3 £ 5 + 

+ 546w 4 t 4 + 510w 5 t 3 + 378m 6 £ 2 + 184w 7 £ + 

+ 41w 8 )M 6 - 2(10t 9 + 55wt 8 + 136w 2 t 7 + 218w 3 t 6 + 265w 4 t 5 + 265w 5 t 4 + 
+ 218m 6 t 3 + 136w 7 t 2 + hhuH + 10w 9 )M 4 + (t 2 + ui + w 2 ) 2 (4t 6 + 22wF + 
+ 37u 2 t 4 + 36u 3 £ 3 + 37u 4 P + 22uH + 4m 6 )M 2 - 



2\4 



+ 2t«(t + U) (f + + u 2 ) 



(14) 



It is clear that, contrary to Xco,2 case, this cross section is finite in small pt region. This is 
explained by the Landau- Yang theorem, that forbids the production of Xci-meson from two 
massless gluons. As a result, the squared matrix element of this reaction is proportional 
to the virtuality of the intermediate gluon (for example, t-channel gluon in figU^i), so this 
factor compensates the divergency, caused by the propagator. 



3.2. gq -> Qq, gq -> Qq 

In this section we will consider charmonium production via scattering of the gluon on 
light (anti)quark. The corresponding diagrams are shown in figure [5l Only Xc mesons can 
be produced in these reactions and the matrix element has the form 

M(gq^xcjq) = ^^Sabuiksh^uih^A^ih-kuh), 

where i,j and a, b are color indexes of quarks and gluons, e v 2 is the polarization vector of the 
initial gluon, u(k 3 ) and u(ki) are the wave functions of quark and antiquark respectively, 
and A$ are the vertices of the interaction of the gluon pair with the Xc-mesons, that were 
introduced earlier. 
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Figure 5: qg — > Qq 



The partonic cross sections of the processes qg — > Xc.il have the form 



di 



d°{qg -> Xc2q) 
di 



di 



8a?7r|i?p 



-^9s 2 M 3 (M 2 -t) 
- 2 t(6s 2 + 7ts + 4t 2 )M 2 + P(2s 2 + 2ts + t 2 



9M 8 - 6(3s + 4t)M 6 + (18s 2 + 30ts + 22t 2 )M 4 



16a 3 7r |i?p 



-^9s 2 M 3 (M 2 -t) 
- 2t(6 s 2 + lis + t 2 )M 2 + P(2s 2 + 2is + i 2 



6M 8 - 12(s + t)M 6 + (12s 2 + 24ts + 7t 2 )M 4 - 



16a 3 7r |i?p 



3s 2 M 3 (M 2 -i)' 
+ i(2s 2 + 2is + t 2 ) 



(4s + t)M 4 - 2(2s 2 + 3ts + P)M 2 + 



We see that, similar to gluon-gluon reactions, that cross section is finite for Q = Xci, and 
diverges in the range of small pt for Q = Xco 2- Their ratio of the divergent cross sections is 



lim 



| da(gq -> Xc2q) / da(gq -> Xeo<?) 



4 

3' 



di I di 

exactly equal to the same ratio in the gluon-gluon case. This fact is not surprising, since 
the nature of the divergencies of gg — > X9 an d gq — > XQ is the same. 



3.3. qq -» Qg 



If we consider charmonia production in pp-annihilation, there are also valence antiquarks 
in the initial hadrons. In this case charmonia production via the quark-antiquark annihi- 



lation is possible. Experimental data says 



231 ] . that the ratio of the -0-production cross 
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sections in pp and pp reactions at the energy of 40 GeV is equal to 

a(pp -*ij) + X) = 
a(pp — > ip + X) 

This is a clear sign of the significance of qq process at low energies. 

This reaction is evidently cross-symmetric with respect to the process considered in the 
previous section, so the matrix elements of these processes can be expressed through each 
other: 



\M(qq^Qg)\ 2 = \M(gq-+Qq)f 



With the help of this relation it is easy to obtain the following expressions for the partonic 
cross sections: 



di 



da(qq -> x*9) 
di 



da(qq -> Xc2g) 
di 



64a 3 yr \R 



■I |2 



27s 3 M 3 (M 2 -s) 4 



9M 8 - 6(4s + 3t)M 6 + (22s 2 + 30ts + 18t 2 )M 4 



2 s(4s 2 + 7is + 6t 2 )M 2 + s 2 (s 2 + 2ts + 2t 2 ) 



128a 3 vr \R' 



[s + At)M 4 - 2(s 2 + 3ts + 2t 2 )M 2 + 



9s 2 M 3 (M 2 -s) 
+ s(s 2 + 2ts + 2?) 
128ct?7r I^R'pl 2 



27s 3 M 3 (M 2 -s) 4 



6M 8 - 12 (s + i)M % + (7s 2 + 24ts + 12t 2 )M 4 



- 2 s(s 2 + 7ts + 6?)M 2 + s 2 (s 2 + 2ts + 2t 2 ) 
These cross sections are finite in whole kinematically available region. 



4. HADRONIC CROSS SECTIONS 

Let us now proceed to experimentally observable values and consider the process 

A( Pl )B(p 2 ) - Q(p)+X, 

where A and B are the initial hadrons, Q = ip or Xcj, and in the parenthesis we show the 
particle momenta. The cross section of this reaction is expressed through the cross sections 
of considered above partonic reactions: 

d<j AB (Q) = da(AB ^Q + X) = J2 d°ab B (Q) = Yl da ( AB ^ ab^ Q + X) = 

a,6,c a, b 

= J dx 1 dx 2 f a /A(xi)fb/B(x2)dcr(ab -> Qc). (15) 
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Here we sum over partons a and b, are the momentum fractions held by these partons, 
and f a /A{ x i)i fb/B(x 2 ) are the distribution functions of the partons in the initial hadrons. 
It is convenient to use other integration variables: 

X = X\ — X2 

and the squared energy of the partonic pair 

s — (h + k 2 ) 2 ~ xix 2 s, 

where s = (pi + P2) 2 and in the last equation we have neglected all masses except the mass 
of the final charmonium. The hadronic cross section is now written as 

(s-M 2 )/(2^i) 

^(Q,p T >A) = [ dp T J g QC) x 



A 



(pr+^P^+M 2 ) 2 



1-s/s 

/dx 
-^fa/A{x l )f b/B (x 2 ), (16) 

-{1-s/s) 



where 



x = x\ + x 2 = \l x 2 H 



s 

the charmonium transverse momentum is defined in eq. (1131) . and 

da 

— + (t^u) 
dt 



da Ptw s 



For px < A it is impossible to distinguish the reactions of the charmonium production with 
the emission of an additional parton and without such an emission. For this reason in small 
Pt region we will consider only LO reactions 

PP -> 99-+ Xc0,2 
and the cross sections of these reactions are equal to 

ab / \ n 2 T(xcj^gg) f dx 

Vlo(XcJ) = — ^3 J -jf g /A{xx)f g /B{x2)- 



In our numerical estimates we used the distribution functions presented in the works 24] 



and [251 ] . Qualitatively these distributions give the same results, though for small energies 
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Figure 6: Cross sections of the charmonium production in proton-proton reaction for different 
values of cm. energy. 1 — total tp production, 2 — tp production through the radiative Xc2 decay, 
3 — direct ifi, 4 — tp produced in radiative Xd decay, 5 — collinear contribution to Xc2 production. 



some quantitative difference is observed. In our paper only the results obtained with the 



help of 



24J are presented. These distributions were calculated at the scale 



Q 2 = 3 GeV 2 



and other parameters are equal to 



a s = 0.3, M = 3 GeV, A = 0.3 GeV' 
\R S (0)\ 2 = 0.81 GeV 3 , 1^(0) | 2 = 0.075 GeV 5 . 



In figure E] we show the dependence of if) production through different processes on total 
energy. The bold line (label 1) shows the summed over all processes 



a = o? g ty) + Br( Xcl -> ^ 7 K P (Xd) + Br( Xc2 - ^ 7 K P (x c2 ), 



dotted line (label 2) — ip production from radiative \c2 decay, direct ip production (thin 
solid line, label 3); V'-production from radiative Xci decay (dashed line, label 3) and collinear 
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Xc2 



Figure 7: The role of the quark-gluon subprocesses versus the energy of the reaction 
contribution (dash-dotted line, label 5). Dots are the experimental results, presented in 



261 ] . It is seen that our calculations give satisfactory agreement with experimental data in 
wide interval of energies, though for small energies there is some difference. This difference 
is caused by the fact, that our cross sections of direct ip and Xci productions are smaller, 
than the experimental results. Another possible reason of the disagreement is that for small 
energies it is not valid to neglect the proton mass in comparison with the energy of the 
reaction. This mass leads to the decrease of the partonic energy and increases the cross 
section and can be taken into account shifting the theoretical curves along the horizontal 
direction. At high energy this effect is not so crucial, as it is clearly seen from figure EJ 

We would like to stress, that the relative contributions of the mechanisms of i/j produc- 
tion depend strongly on the energy of the reaction. For example, in figure [7] we show the 
dependence of the ratio 



Br( Xcl -> ^ 7 )[<(Xci) + gg(x c i)] + Br( Xe2 -> ^)[af u ( Xc2 ) + <g(x c2 )] 



on the energy. This ratio describes the role of quark-gluon subprocesses in the production 
of -0-meson. In is seen, that it depends strongly from the energy, but its value is suppressed 
in the whole shown range. 

Our results for Xci and direct ip production cross sections are somewhat smaller, than 
the experimental values. From experimental data [3( we know, that the ratio of Xci and Xc2 
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Figure 8: The ratio of the cross sections of ip production in pp and pp reactions 
production cross sections is equal to 

4^ - 0.3, 

while our calculations give the value approximately 3 times smaller. We think that the 
source of this difference is that only one of the producing Xci gluons was virtual. According 
to figure [2] the axial charmonium is very sensitive to the virtualities of the gluons, so this 
effect will increase the cross section of its production. The Xco an d Xc2 production cross 



27] it was shown, that the 



section would not change significantly. As for the direct ip, in 
process qg — > ipqg almost doubles the value of its production cross section. It should be 
noticed, however, that the collinear singularity, appearing in this process, was removed by 
introducing the mass of the quark, and the value of this mass was chosen to be m q = 5 
MeV. We think, that this value is too small, and hence, the role of qg — > ipqg process is 
overestimated in 271 ]. 



In the case of pp annihilation there are also valence antiquarks in the initial hadrons. 
That is why the cross section of pp — > Q + X reaction will receive the contributions from the 
quark- ant iquark subprocesses, that were considered in section 3.3. In figure [8] we show the 
dependence on the energy of the ratio of the cross sections a(pp — > ipX) and a(pp — > ipX). 
It can be seen, that this ratio is significant for small energies and tends to unity for higher 
ones. The reason for this is that the quark-antiquark cross sections are suppressed by the 
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Figure 9: Transverse momentum distribution of the total -0-production cross section at energy 
yfs « 27 GeV. Dashed line — our cross section, solid line — smeared distribution, dots — experi- 



mental values 



28(]. 



propagator of the s-channel gluon. Such behavior of this ratio is observed in the experiment. 
Staring from formula ( fT5l) it is easy to obtain the distribution of the charmonia production 



now 



28] 



cross section over the transverse momentum. In figure [9] with the dashed line we s 
this distribution for energy -y/i = 27 GeV and the experimental data, taken from 
It is clear, that our technique gives too fast decrease of the distribution with the rise of 
transverse momentum. In order to remove this disagreement one can take into account the 
transverse motion of the initial gluons in the proton. This method is used in the so called 
factorization to explain the high energy results. For low energies, unfortunately, it cannot 
be used, since the unintegrated distribution functions at these scales are poorly known. We 
will take into account the transverse motion of initial gluons by phenomenological correction 
[2 



da A T , da 



dpr Ana 2 J dq> 

where 



A = 3, a = 0.5 GeV. 
The corrected distribution is shown in figure [9] with solid line. 



(18) 
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1 da nb 
px dpj- ' GeV 

10 p 




1.25 1.5 1.75 



p T , GeV 



Figure 10: Transverse momentum distributions of tp production. Thick solid line — total cross 
sections; thin solid line — direct ip production; dashed line — ip form radiative Xci decay; dotted 
line — ip form radiative \c2 decay 

4.1. £ = 70 GeV 



In this section we will consider in more details charmonium production at 8 = 70 GeV. 

In figure [10] the Pt distributions of ip productions form different mechanisms are shown: 
direct ip production (thin solid line, label 3), ip form radiative Xc decays (xci — dashed line, 
label 4; \c2 — dotted line, label 2) and total cross section (thick solid line, label 1). We 
have not performed any smearing, since the value of the parameters a and A is unknown. 

The other widely used kinematical variable is the longitudinal charmonium momentum 
Pl- It is, however, more convenient to work with dimensionless variable 

2p L 



that can be expressed through the introduced earlier values: 

r2 s M 2 _ 



Xp 



2s 



M z s 
x + 



where 



cos 6 
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1 + 



2s 



2t 



-xz, 



(19) 



s-M 2 ' 
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dcr 

dxp ' 
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-0.2 



0.2 
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Figure 11: xp distributions. Solid bold line — total cross section; thin solid line — direct ip 
production; dashed and dotted lines — V from Xci and Xc2 respectively 

and #i3 is the angle between the momenta ki and k3, measured in partonic rest frame. In 
figure [11] we show x^-distributions of the cross sections for ip and x states. It can be seen, 
that radiative Xc2 decay gives the main contribution to ip production cross section. 

For the determination of the gluon distribution functions from the experimental data one 
needs to know the momentum fractions of these partons. These fractions, unfortunately, 
cannot be determined directly. So, it is necessary to know what values of x\$ correspond to 
specific value of the variable xf- If the energy of the partonic reactions coincides precisely 
with charmonium mass (that is we deal with the resonant charmonium production), than 
only the first term of relation (fT9l) gives nonzero contribution and we have one to one 
correspondence between x li2 and xp'- 



If the energy of the partonic reaction is larger than the charmonium mass, this correspon- 
dence is ambiguous. We know, however, that the states with large partonic energies are 
suppressed by the partonic luminocities (that is the probability to observe partons with 
large momentum fractions x\ and x^), so the contribution of these states into cross sections 
will not be significant. For this reason the errors caused by the mentioned ambiguous will 
be small. In figure [12] we show the correspondence between x\p and xf and the errors in 




(20) 
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Figure 12: Correspondence between experimentally observed variable Xf and momentum fractions 

^1,2 

the determination of these variables. The dashed lines correspond to the relation fTSOj) with 
s = 1.3M 2 . The x% error was determined from the "half-width" condition, that is 

d 2 a 1 (fa 

r 1 max 

(the same for x-z). We would like to mention, that these errors have nothing to do with 
the experimental uncertainties and cannot be reduced by the increase of the experimental 
precision. 



5. SPIN ASYMMETRIES 



The ultimate goal of the PAS-CHARM experiment is the determination of the spin dis- 
tribution functions of the partons in the polarized nuclon f(X,x). Here A = +1 corresponds 
to the case of the parton spin aligned along the spin of the parent nucleon, while A = — 1 
corresponds to opposite alignment of these spins. For the determination of these functions 
it is proposed to consider the spin asymmetry 

A _ <r(tT)-<r(U) (21) 

where symbols j and j in the ^-production in the proton-proton reaction cross section 
correspond to the cases with proton polarized along the direction of its motion or against 
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this direction. 

The explicit expression for asymmetry (|2T]) is 

_ A A/(xi) A/(s 2 ) 

where Af(x) = f(+,x) — f(—,x) and the partonic spin asymmetry is connected with the 
cross sections of charmonium production from two partons with the Helicities Ai and A 2 : 

da(++) - da(+-) 



da(++) + da(- 



Using the helicity cross sections of charmonium production in gluon-gluon annihilation 22 ] 
and the experimental constraints on the ratio 

AG = Af g/P (x) 

G fg/ P (x) 
[ — I 

that reaches its maximum value 0.9 for x ~ 0.3 [30(] it is easy to obtain the estimate for 
asymmetry All a t energy 70 GeV. In figure [13] we show the pt distributions of the partonic 
asymmetries (111- It is clearly seen, that different charmonia states give the opposite sign 
contributions into this asymmetry, so averaging over these states leads to significant decrease 
of the asymmetry. So it is necessary to determine the type of the experimentally observed 
charmonium meson and separate the contributions from directly produced tp and ip from 
radiative Xci,2 decays. 



6. CONCLUSION 



In the proposed experiment SPASCHARM it is planned to measure the polarization 
distribution functions of quarks and gluons in the nucleon. For this the detailed analysis 
of different modes of charmonia production in the hadronic experiments is necessary. The 
collinear approximation, that is widely used in parton model, suffers from a number of serious 
drawbacks, that make it impossible to use this approximation in our case. First of all, this 
approach gives no information on the distributions over the transverse momentum of the 
final charmonium. These distributions give rich information on the charmonia production 
mechanisms. The second drawback is that at leading order in the strong coupling constant 
only x c0 2 -mesons can be produced, while the experimental data gives large direct if) and Xci 
signals. 
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Figure 13: px distributions of the partonic spin asymmetry o,ll f° r reactions gg — > (solid line), 
55 — * Xci5 (dashed line) and gg — > x C 2<7 (dotted line) 

The natural way to solve these difficulties is to consider the third-order processes in 
the strong coupling constant. This is the topic of our article. We have shown that the 
consideration of the higher order processes increases significantly the number of mechanisms 
of charmonia production. In addition to Xco,2 production in gluon-gluon fusion, that is 
possible in the leading order, there appears also a production of direct ip and Xci an d 
production of P-wave charmonia in quark-gluon interaction gq — > XcjQ- in the case of proton- 
antiproton annihilation the production of x c - m esons in the process qq — > Qg is possible, that 
explains the experimentally observed excess of charmonia production in proton-antiproton 
annihilation at low energies in comparison with the proton-proton scattering. 

The results, presented in our paper show, that this approach gives the cross section of 
total ip production in proton-proton interaction, that describes well the experimental data 
in a large range of energies, and qualitatively describes the xa production and distributions 
over the transverse momentum of the final charmonia. There is, however, no quantitative 
coincidence. Our values for Xci production cross section are about 3 times smaller, than 
the experimental values, and pt distributions fall with the transverse momentum faster, 
than it is observed experimentally. It is known from the theoretical and experimental anal- 
ysis of TEVATRON data, that these contradictions can be removed with the help of kr- 
factorization. Unfortunately, For low energies (for example the energies of U70) it is hard 
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to use this approach, so we have introduced phenomenological corrections. 

The presented in our paper detailed analysis of different mechanisms of charmonia pro- 
duction in proton-proton interaction at the energy of the initial proton equal to 70 GeV 
shows, that the relative contributions of these mechanisms depend strongly on the trans- 
verse momentum of the final charmonium and other kinematical variables. 

In the last section of our article we have given the expression for spin asymmetries of 
partonic processes of charmonium production in gluon-gluon fusion. These asymmetries, 
with the analogous values for quark-gluon subprocesses, are necessary for the determination 
of the polarization distribution functions of quarks and gluons in proton. From our results 
it is clearly seen, that the contributions of different charmonia states into this asymmetry 
compensate each other. For this reason it is necessary do separate directly produced if) 
mesons and ?/>-Hiesons produced via radiative Xci,2-decays. 

The authors would like to thank A.N. Vasiliev and V.V. Braguta for useful discussions. 
Our work was supported in part by Russian Foundation of Basic Research under grant 
07-20-0041 7a. 
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